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ABSTRACT 


The fundamental mechanics of delamination in fiber composite laminates is 
studied. Mathematical formulation of the problem is based on recently developed 
laminate anisotropic elasticity theory and interlaminar fracture mechanics con- 
cepts. Stress singularities and complete solution structures associated with 
general composite delaminations are determined. For a fully open delamination with 
traction-free surfaces, oscillatory stress singularities always appear, leading to 
physically inadmissible field solutions. A refined model is introduced by con- 
sidering a partially closed delamination with crack surfaces in finite-length 
contact. Stress singularities associated with a partially closed delamination 
having frictional crack-surface contact are determined, and are found to be dif- 
ferent from the inverse square-root one of the frictionless-contact case. In the 
case of a delamination with very small area of crack closure, a simplified model 
having a square-root stress singularity is employed by taking the limit of the 
partially closed delamination. The possible presence of logarithmic— type stress 
singularity is examined; no logarithmic singularity of any kind is found in the 
composite delamination problem. Numerical examples of dominant stress singular- 
ities are shown for delaminations having crack-tip closure with different fric- 
tional coefficients between general 0 ^ and 02 graphite-epoxy composites. 
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1. INTRODUCTION 


Delamination has been a problem of significant concern in the reliable 
design and analysis of advanced fiber composite laminates . Separation of 
composite laminae, caused by high local interlaminar stress and low strength 
along the ply interface, can result in destruction of load transfer, reduction 
of stiffness, and loss of structural integrity, leading to final structural 
and functional failure. From the mechanics point of view, delamlnatlon 
involves initiation and growth of macroscopic cracks between dissimilar, 
strongly anisotropic solids. A rigorous mathematical study of delamlnatlon is 
recognized to be difficult, especially in a finite-dimensional fiber composite 
laminate. The complexities include the inherent crack-tip singularity, the 
effect of anisotropy of each constituent fiber lamina, and the abrupt change 
of stiffness or ply orientation through the laminate thickness direction. In 
addition, the three-dimensional state of stress and deformation associated 
with the composite delamination always gives rise to a combined opening (mode 
I), in-plane shearing (mode II), and out-of-plane tearing mode (mode III) 
fracture, which render the problem mathematically intractable in many 
cases. The mechanics of delamination in fiber composite laminates is, 
therefore, not only of significant academic interest but of practical 
importance. In this paper, the first of two articles in a row, the 
fundamental nature of stress singularities and associated field solutions for 
a delamination in a fiber composite laminate are investigated. 

Owing to the aforementioned complexities , studies on an Interface crack 
between dissimilar anisotropic materials have been limited. Gotoh [1] appears 
to be the first to examine the two-dimensional problem of partial debonding 
between dissimilar anisotropic plates under a plane stress condition. 

Clements [2] has used Stroh's approach [3] to study the problem of an 
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interface crack between two generally anisotropic half -spaces. Willis [4] has 
also conducted a two-dimensional stress analysis of a crack on the plane 
interface of two bonded dissimilar half -spaces. The analysis has been 
combined with the usual local form of Griffith's virtual work argument to give 
a failure criterion, involving a stress concentration vector and specific 
surface energy of the bonded interface. All of the asymptotic solutions 
obtained in [1,2,4] have an oscillatory displacement field that material 
interpenetration on either side of the crack surface is predicted. Similar to 
those found in the solutions for an Interface crack between dissimilar 
isotropic materials [5-9], these physically unreasonable results have led to 
the argument of solution inadmissibility for the crack problem in dissimilar 
anisotropic media. To correct the unsatisfactory feature of oscillatory 
stress singularity, Wang and Choi [10] have recently reconsidered the problem 
of an interface crack between dissimilar, strongly anisotropic fiber-composite 
^^^^~spaces by introducing a partially closed interface crack model, in which 
the crack is not completely open and that its surfaces are in frictionless 
contact near the tip. The formulation leads to a singular integral equation, 
which is solved numerically. Numerical results from this refined model [10] 
exhibit an inverse square-root stress singularity and, therefore, physically 
meaningful fracture mechanics parameters can be defined consistently with 
those in fracture problems of homogeneous materials [11-13] and in the model 
given by Comninou [14,15] for an interface crack between two isotropic 
media. Moreover, significantly global crack closure has been found [16] for 
an interface crack between dissimilar anisotropic elastic half-spaces 
subjected to mixed-mode loading — a situation that is generally experienced by 
a delamination in finite dimensional fiber composite laminates. 
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In this paper, we employ Lekhnltskii *s complex-variable stress potentials 
[17] in conjunction with an eigenfunction expansion method to examine the 
mechanics and the mathematical solution structure for a delamlnatlon with 
frictional crack-tip closure in a composite laminate. Based on the general 
solution structure determined, an advanced numerical method using singular 
finite elements is then developed to study the detailed delamlnatlon behavior 
in finite dimensional fiber composite laminates with any arbitrary 
combinations of lamination parameters, geometric variables, and crack 
dimensions. Owing to space limitation, the numerical method and the detailed 
composite delamination behavior are reported in an accompanying article [18]. 

In the next section, the problem definition and basic assumptions are 
stated. Basic laminate anisotropic elasticity equations and formulation of 
the composite delamlnatlon problem are Introduced in Section 3. General 
solution structures for asymptotic stress and displacement fields are 
obtained. Stress singularities associated with an open interlaminar crack and 
with a partially closed delamination tip with frictional crack-surface contact 
are determined respectively in Section 4. Influences of frictional 
coefficients on delamination stress singularities are examined. A simplified 
model for a delamlnatlon with a very small area of crack-tip closure is also 
introduced. The possibility of existence of additional singularities in 
logarithmic forms in homogeneous and particular solutions is investigated. 
Results are presented for delaminations with different local crack-surface 
traction boundary conditions in composite laminates containing various fiber 
orientations. The eigenvalues and associated stress singularities obtained in 
this study provide the most fundamental information on complete solution 
structures of delamlnatlon stress and deformation fields, and establish a 
basis for formulation of the singular finite elements used in the next paper 
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[18] to study the detailed delamlnatlon behavior in finite dimensional 
composites with general lamination and geometric variables. 
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2. STATEMENT OF THE FROBUM AND ASSUMPTIONS 

The problem considered here is a composite laminate (Fig, 1) composed of 
unidirectional fiber-reinforced plies of uniform thicknesses, hj, h 2 » h 3 , 
hjj. The composite has a finite dimension with a width equal to 2b. For 
simplicity and without loss of generality, we restrict our attention to the 
cases of symmetric composite laminates with fiber orientations 
[0 j/02/93/..»03/02/®i 1* thicknesses are also symmetric with respect to 

the x-z plane, l.e,, for each ply above the x-z plane (y > 0), there exists an 
identical ply with the same ply thickness below the x-z plane (y < 0), 
Delamination with a length a is assumed to occur in the form of an Interface 
crack between dissimilar, strongly anisotropic fiber-reinforced composite 
laminae with fiber orientations 0jj and 

The composite laminate is assumed to be subjected to tractions acting in 
planes normal to the z-axis and distributed uniformly along the z-axis without 
variation. In the case that the finite dimensional composite laminate has a 
finite length, axial loads and moments are assumed to act on the ends of the 
composite body. The composite laminate is further assumed to be sufficiently 
long that in the region away from the ends, end effects are negligible by 
virtue of the Saint Venant principle. Consequently, the components of 
stresses in the laminate are independent of the z-axls. The special case in 
which all components of stresses and displacements in the composite are 
Independent of the z is well-known as the generalized plane deformation 
problem [17], 

The objectives of this study are to: (1) establish a mathematical basis 

for the mechanics of delamination based on laminate elasticity theories and 
Interlaminar fracture mechanics concepts; (2) determine stress singularities 
and associated solution structures for composite delaminations with different 
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local crack-tip deformation configurations; (3) obtain asymptotic stress and 
deformation fields governing the fundamental behavior of delamination; and 
(4) study the Influences of various lamination and material variables such as 
fiber orientation and crack-surface frictional coefficients on the 
delamlnatlon stress singularities. 
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3. LAMINATE ELASTICITY BASIS FOR (XRfPOSITE DELAMIHATION 

The development of the mechanics of composite delamination is based on 
recently established theories of anisotropic laminate elasticity [19,20,21] 
and fracture mechanics concepts of interface cracks between dissimilar, 
strongly anisotropic composites [4, 10, 16,21] • In this section, governing 
partial differential equations for composite laminate elasticity problems are 
established first, based on Lekhnitskii *s complex-variable stress potentials 
[17]. General solutions for the laminate elasticity problem with interlaminar 
cracks are introduced. Stress singularities associated with a composite 
delamination having homogeneous local boundary conditions are defined. 

Solution structures of asymptotic stress and displacement fields are 
constructed for a delamination between dissimilar general fiber composite 
laminae. Additional terms of logarithmic forms in the homogeneous and 
particular solutions for the composite delamination problem are examined. 

3.1 Basic Equations 

The fundamental mechanics of delamination in a fiber-reinforced composite 
laminate may be studied from the schematics illustrated in Fig. 1. The 
constitutive equations of each fiber-reinforced composite lamina with 
rectilinear anisotropy of a general form in the structural (x-y-z) coordinates 
are denoted by generalized Hookers law in contracted notation as 

Oj (i,j=*l,2,3,...,6), (1) 

where the repeated subscript indicates summation, and is a compliance 
tensor. The engineering strains in Eq. 1 are defined in a Cartesian 
coordinate system by 
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where u, v, and w are displacement components. The stresses are defined in 
an analogous manner in the Cartesian coordinate system# For a composite 
laminate in the aforementioned loading condition, mathematical formulation for 
this class of elastostatic problems can be made using the well-known 
Lekhnitskii complex-variable stress potentials [17], F(x,y) and ^(x,y), 
defined as 

3^F _ _ 32p 
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Following the same procedure in [17,19], we can easily obtain the following 
system of partial differential equations for each anisotropic composite 
lamina : 


L 3 F + l.^'V - 2 A^ + - A 2 S 33 , 

(4a) 

L^F + = 0 , 

(4b) 


where L£, Lg and are differential operators of the second, third, anH 

fourth orders which have the form: 
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with 


\j ■ ^i3 ®3j^^33 


(i.j = 1,2, 4, 5, 6), 


( 6 ) 


The constants Aj and A 2 in Eqs. 4(a) and 4(b) characterize the bending of the 
composite body in the x-z and y-z planes respectively, and is the relative 
angle of rotation about the z-axis. 

Assume that external tractions on the lateral surface of the composite 

cross section are given as T , T and T . The boundary conditions on the 

X y z 

lateral surface bB are as follows: 

a n + T n = T , 

XX xy y X 


^xy ^x ^y ^y " *^y* 


(7) 


^xz *^x ^ ^yz ^y *^z* 


where n^ are directional cosines of the bounding surface 5B. The conditions 
at the ends of the composite have the form: 


B ^ B ^ B 

// y dx dy = M^, // x dx dy = M 2 , (8) 

B B 

// (Ty2 X - 
B 

where the integrals are taken over the entire area of the cross section B, and 
Mj, M 2 , and M^ are applied force, bending moments, and twisting moment at 
ends of the composite, respectively. 


3.2 General Solutions 

The general solutions for the governing differential equations have been 
shown [17] to have the form as 
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6 

■' - Jj (5) 

6 

( 10 ) 

where the complex variables are defined as Zj^ = x + y; and are 
particular solutions of the nonhomogeneous system; the prime (*) denotes 
differentiation of the analytical functions Fj^ with respect to their 
arguments, and are roots of the following algebraic characteristic 
equation : 



-^4(fi)A2(p) - A§(p) = 0, 

( 11 ) 

with 

^ 2 ( 1 ^) “ ^ 55 ^^^ “ 2 S^^p + 

( 12 a) 


^3(1.) - -S,5i.3 - (§3^ + S5^)^2 + (§25 + 3^5)^ - ■§2^. 

( 12 b) 


.,(k) = 3„." - 23^5,3 ^ ^ ^ - 3 ^^. 

( 12 c) 


The Tij^ in Eq. 10 are complex numbers equal to 


( 13 ) 

We now choose the form of Fj^(Zj^) as 

= Sc Z^^/[(S+1)(5+2)], ( 14 ) 

where and 6 are arbitrary complex constants to be determined. Substituting 
Eq, 14 into Eqs. 3, 9 and 10, we obtain the homogeneous solutions for stress 
and displacement components in polar coordinates (r-<l)-z) as follows: 
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where 


= r(e^'t’ + \^e ^‘*’)/(l + 


= (1 + iyj^)/(l - i\)* 


(17a) 

(17b) 


The coefficients Hj^]j^('}>) (i = 1»2,3 8; k = 1,2,3) are known functions 

of ^9 and , defined in Appendix 1. 

The complete laminate elasticity solutions for the composite mechanics 
problem can be written as 


“3 - + “j”' . 


(i = 1,2, 3, 4, 5, 6), 
(j = 1.2,3), 


(18a) 

(18b) 


where and are particular solutions associated with the loading 

condition of each individual case studied. The expressions for and 

z z 


can be obtained as 
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and 


„Ch) . 


S3J o(«/S33, 



(AjX + A^y + A^) 




/S 


33 


(18c) 


(J “ 1»2»4,5,6). (18d) 


3»3 Asymptotic Stress and Displacement Fields 

Using Eqs. 15 and 16 and applying local homogeneous traction boundary 
conditions on crack surfaces 58^. (Fig. 2) and interface continuity (matching) 
conditions along we obtain a system of twelve homogeneous linear 

equations in =m, m + 1), i.e,, 

D C = 0, (19) 

where D(6) is a 12 x 12 coefficient matrix involving 6 in a transcendental 
form, and C is an unknown 12 x 1 column eigenvector. The nontrivial solution 
U requires that the determinant of the coefficient matrix vanishes, i.e., 

IID(6)II = 0. (20) 

This leads to a standard eigenvalue problem, and the 6 can be determined from 
the transcendental characteristic equation. Standard numerical methods such 
as the Muller method [22] with the aid of a digital computer are needed for 
this purpose. The eigenvalues determined from Eq. 20 provide important 
information on the fundamental structure of stress and displacement solutions 
for the composite delamlnatlon problem. Furthermore, the eigenvalues 6^ which 
satisfy the following condition: 

-1 < Re [6^] < 0 (21) 

characterize the fundamental nature of stress singularities and provide the 
asymptotic stress and deformation fields at the delamination tip. In the case 
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that local crack surface tractions are nonvanishing, for example, the crack 
closure problem, Eq. 19 needs to be modified. Delaminations having crack-tip 
deformation configurations with nonvanishing local traction boundary 
conditions are discussed in detail in the next section. 

For a delamination problem in composite laminates with general lamination 
variables and fiber orientations, the algebraic multiplicity of the 
eigenvalues determined from Eq. 20 may give rise to additional terms of the 

6n 

logarithmic form (An Z^) in the homogeneous solution, as first suggested 
by Dempsey and Sinclair [23,24]. In this situation, the following terms may 
also be a part of the homogeneous solution in addition to Eqs. 15 and 16: 

6 ^m 6 

^ ^ik “ 1,2, 4, 5, 6), (22a) 

56 “ ■ 
n 

6 J^m (6 + 1 ) 

"J - I \ 

56 “ 
n 

where is the order of the logarithmic multiplier in the eigenfunction 
corresponding to eigenvalue 6^, and is related to the property of the D 
matrix hy - M-(N-R), in which M is the algebraic multiplicity of the root 
5^, and N and R are the order and the rank of the D matrix, respectively. The 
presence of the logarithmic terms, Eqs. 22(a) and 22(b), in the homogeneous 
solution requires a nontrivial solution for Cj^. Detailed discussion of the 
conditions for the existence of Eqs. 22(a) and 22(b) in the composite 
delamination problem can be found in [25]. 

In the construction of asymptotic solutions for delamination stress and 
displacement fields, the particular solution for the system of governing 
differential equations also contributes to the complete solution. It is 
apparent that the structure of the particular solution is related to the 
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applied loading and deformation of the delaminated composite. For the 
convenience of further developments, we consider here the case of a composite 
laminate with delamlnatlons subjected to a uniform axial stretching, l.e., 

“ S* Under this circumstance. It has been shown [25] that the particular 
solution has a similar form as Eqs. 22(a) and 22(b), 


,(P) = 


“oi - 

d6 ° " 

n 


(1 = 1,2, 4, 5, 6), (23a) 


“ “oj =0 = 1*2,3), (23b) 

56 ° n 


-9 ° V „(p) 


where the components and u^j are known quantities determined by the remote 
loading condition. The In Eqs. 23(a) and 23(b) Is the order of the 
logarithmic eigenfunction at 6^^ = 0 and Is related to the multiplicity of 
the root = 0 and the rank and the order of the matrix D by = Mq-(N-R). 
[Note that Eqs. 23(a) and 23(b) contain logarithmic terms of the forms 


(In Zj^), (JLa Zj^)2, ... (Jtn Zj^) °.] The necessary and sufficient conditions 
for the existence of the particular solution, Eqs. 23(a) and 23(b), can be 
shown [25] as 


r*(L) * 

^h P “0 (24) 

for every left eigenvector of D*(0) defined In [25], where p* Is a 

loading vector resulting from and u^ j , and the dot (•) denotes the inner 

product of the two column vectors. In the case that Eq. 24 does not hold, one 

needs to consider the logarithmic terms of a higher order through a hlgher- 

+1) (A +1) 

order differentiation 9 /56^ ° in Eqs. 23(a) and 23(b). A detailed 

discussion on this Is given in Reference 25. 
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4« DELAMIMTION STRESS SINGULARITIES IN COMPOSITE LAMINATES 


Based on the general solution structures given in Eqs. 15, 16, 22 and 23, 
it is possible at this point to examine the detailed nature of stress 
singularity associated with a delamination in a fiber-reinforced composite 
laminate# Because of the local nature of the stress singularity, we focus our 
attention on the crack-tip region of a delamination between the mth and 
(m+l)th laminae with fiber orientations 9^ and respectively# Both fully 

open and partially closed delaminations are considered# In the case of a 
delamination with an extremely small area of crack-tip closure, a simplified 
model by taking the limiting case of a partially closed crack is introduced. 


4.1 Delamination with Traction-Free (Fully Open) Surfaces 

Assuming that the crack surfaces are fully open and that the interface 
5 B]q between the plies is perfectly bonded along r > 0 as shown in Fig# 2, we 
can immediately introduce the local traction-free boundary conditions along 
the delamination surfaces <t) = ± 11 , 






(25a) 

(25b) 


The continuity (or matching) conditions of interlaminar stresses and 
displacements along the ply interface (J) = 0 read as follows: 




(m+l), 


.(m+1). 




More explicitly, the local homogeneous boundary conditions, Eqs, 25(a) and 
25(b), and continuity conditions, Eqs. 25(c) and 25(d), have the forms as 
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I {C<"> [o'-'ot)]®} . 0, 

k=l 


( 26 a) 


k=l 


.(m+1), xi6 , p(m+l) „(m+l) , 




J^(«4">r<f . c« T<f - o. 

(i = 1,2,3; j = 1,2, 3, 4, 5, 6), (26c) 

where 

e~^‘'’)/(l + \^“^), (26d) 

and 


Ta(a)_ 

Ik 


1 , 



r 


(“)= ,/“) 
3k ^ ’ 




T.(a)_ 

6k " 


t 


(a) 

k 


(a = m, m+1). (26e) 


Equations 26(a), 26(b) and 26(c) consist of a system of twelve 
homogeneous linear algebraic equations in The existence of a 

nontrivial solution requires that Eq. 20 holds, leading to a standard 
eigenvalue problem. The solution for 6^^ can be obtained easily from Eq. 20 
and shown to have the form, 

= (n - V 2 ) ± iy> (n ~ V 2 ) , and n, (27) 

where n = 0,1,2,...“>; and y is a constant related to material elastic 
properties, of the adjacent mth and (m+l)th fiber-composite 

laminae. In general, the value of y needs to be determined numerically from 
Eq. 20, which involves in a transcendental form. It is Important to note 
that the eigenvalues of 6^^ obtained from Eq. 27 give critically Important 
information on the fundamental structure of stress and displacement solutions 
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for delaminated composite materials. We remark that the eigenvalues 

(n- V 2 ) i iy arid ("" V 2 ) ate single roots and that all the integers, n, 

including zero, have an algebraic multiplicity of 3 in general. As mentioned 

in the preceding section, the values of which meet the condition Eq, 21 

provide the exact strength (or order) of the inherently stress singularity for 

the asymptotic stress solution at the delamlnatlon crack tip. The possible 

presence of weaker singularities and related terms in logarithmic forms as 

JL 6 i 

discussed in Section 3,3, i.e., (in Zj^) ° and/or Zj^'^(Jln Zj^) in the 
homogeneous solution as well as in the particular solution will be discussed 
later. 


For the purpose of illustration, consider a delamlnatlon located along 
the ply interface of 9/90° graphite-epoxy composite laminae (Fig, 2), The 
following material elastic constants of high-modulus unidirectional graphite- 
epoxy are used in the computation: 

El = 20 X 10 ® psi (137.9 GPa), 

Ej = E 2 = 2.1 X 10® psi (14.48 GPa), 

Glt = Gtz = Glz = 0.85 x 10® psi (5.86 GPa), (28) 

''LT “ ''Lz ~ ''Tz ” 0.21, 


These ply elastic constants are used in the computation here only to 
illustrate the general nature of the current problem* (These constants are 
selected for historical reasons because they have been used in many previous 
studies of the mechanics of composite laminates [19 >26, 27]*) Numerical 
results based on real material constants of the commonly employed T300/5208 
graphite-epoxy with 

El = 19.5 X 10® psi (134.45 GPa), E^ = E^ = 1.48 x 10® psi (10.2 GPa), 

*^T = ®Lz = 0*8 X 10® psi (5.52 GPa), = 0.49 x 10® psi (3.38 GPa), 

^LT = ^Lz = 0.3, = 0.49, (28a) 

are also given in Tables 1 and 2 for comparison. The differences between the 
two cases are generally very small. 
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where the subscripts L, T and z denote the fiber, transverse and thickness 
directions of the composite lamina, respectively. The first three eigenvalues 
5i(i = 1»2,3) which satisfy the aforementioned constraint condition Eq, 21 are 
given in Table 1 to illustrate the exact strength of the stress singularity 
associated with the delamination. To demonstrate further the general 
characteristics of the stress singularities for delamination, results for an 
interlaminar crack between 30°/9 graphite-epoxy composites with the same ply 
properties are also shown (Table 2) for various fiber orientations 9’s. From 
Tables 1 and 2, we observe that a fully opened delamination between dissimilar 
highly anisotropic laminae always possesses three distinct stress 
singularities, i.e., a pair of complex conjugates, 2 “ iy » and a 

real constant, 63 = -0.5. This situation is unique and apparently different 
from the cases of an interface crack between two dissimilar isotropic media or 
orthotropic solids in that the three distinct dominant stress singularities, 
6 j, 62 and 63 , always exist simultaneously in the present fiber-composite 
delamlnatlon problem. In the special cases when a delamination is located in 
the 90“/90® or 30'’/30“ composite system the classical inverse square-root 
singularity for crack-tip stresses is fully recovered as shown in the Tables, 
because the composite laminate becomes unidirectional. We note here that the 
imaginary part of 6 j and 62 , i.e., the value of y# is generally very small as 
compared with the real part of in all cases studied. 

4.2 Delamlnatlon with Crack-Tip Closure 

From Eq. 27 and from the results shown in Tables 1 and 2, it is clearly 
seen that the asymptotic delamination stress field in dissimilar anisotropic 
composites possesses the well-known oscillatory singularities. The associated 
displacement field also exhibits oscillatory characteristics, leading to 
controversial crack-surface interpenetration or overlapping, which is 
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physically inadmissible. Similar results have also been noted by several 
investigators in studying an interface crack between dissimilar isotropic 
materials. In recent studies, Wang and Choi [10,16] have shown that for a 
delamlnatlon between dissimilar, strongly anisotropic fiber composites with 
certain combinations of laminar elastic properties, ply orientations, and 
loading conditions, global crack surface closure may occur. Under these 
circumstances. Interlaminar crack-surface contact or closure needs to be 
considered. 

Consider the case that a delamination a is located between the mth and 
(m+l)th laminae and a portion of the crack surface, c, is closed as shown in 
Fig. 3. Frictional coefficients associated with interface 

([) = 0 are denoted by and f^.^, respectively. An exact analytical complete 
elasticity solution for the delamination problem with crack closure is 
generally difficult to obtain because the unknown contact stress distributions 
along the crack-closure region need to be determined as a part of the final 
solution [10,14], However, the local stress singularities, asymptotic field 
solutions, and associated characteristics can still be determined exactly by 
using the same procedure discussed in Section 4.1 but with some 
modifications. Referring to Fig. 3 for a partially closed delamination. 
Instead of using Eqs. 25(a) and 25(b) we can Introduce the local boundary 
conditions in the crack-closure region (- c < r < 0) as follows: 


oS^(r,TT) - a'[®'^^\r,-Tr) = 0, 






u^“^(r,TT) - u^“"^^\r,-ir) = 0, 
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Along the ply interface (|) = 0, the same continuity conditions Eqs. 25(c) and 
25(d) for and (a = m, m+1) are applicable. 

Using Eqs. 25 and 29 and following the same procedure given in Section 
4.1, we can immediately determine the eigenvalues for a composite 
delamination with crack-tip closure. The numerical example of a delamination 
located along the interface of 0/— 0 graphite— epoxy composites is studied here 
first. Stress singularities associated with the partially closed delamination 
crack tip with different values of frictional coefficients and are 
shown in Table 3. The crack-tip stress singularity is found to be always -0.5 
with an algebraic multiplicity of 2 (i.e., double roots = 62 = -0.5) for 
the delamination having crack surfaces in frictionless contact (i.e., = 

/(j)Z - fact, the inverse square-root stress singularity, = 62 = 

—0.5, is found for all delaminated 0/— 0 fiber composites with frictionless 
crack-surface contact. In the cases of * 0 and/or * 0, stress 
singularities always possess an invariant constant 6^ = -0.5 (single root) as 
in the aforementioned frictionless contact case, and a 62 (with 62 * -0.5, 
single root), which depends on values of the frictional coefficients (Tables 3 
and 4). In Table 3, values of 62 for all delaminated 0/-0 fiber composites 
studied are observed to be slightly larger than -0.5, when > 0 is 
considered. That is, frictional contributions lead to a weaker delamination 
stress singularity 62 than that in a frictionless contact case and in a 
conventional homogeneous open crack case. We note here that for a 
delamination between 0/-0 graphite-epoxy composites, the stress singularities 
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are always Independent of the value of the frictional coefficient because 

of the s)nmnetry (and antisymmetry) of components in elastic stiffness matrices 
of the 0 and —0 plies and the decoupling of from Oy and Xy^ in the 
formulation. This phenomenon is clearly seen in Table 3, where dj is always 
-0.5 and 62 differs from -0.5 gradually as the value of increases. Note 
further that deviations of 62 frictional contact cases from the 

conventional square-root singularity are rather small for all 0/-0 graphite- 
epoxy delamlnatlon problems. 

Stress singularities are also determined for delamlnatlons in more 
general cases of 0^/02 fiber composites with 0^ ^ 02 * For illustration, the 
results of a delamination between 30°/90® graphite-epoxy composites are 
presented in Table 4 for various values of and seen from the 

Table that the double roots = 62 = -0.5 also appear for a delamlnatlon in 
0^/02 composites with crack surfaces in frictionless contact 
However, in the cases of a delamination with crack surfaces in frictional 
contact, 62 is apparently influenced by the values of both and In 

Table 4, values of 62 for the partially closed delamination in 30‘’/90'* 
graphite-epoxy composites with different and observed to be 

smaller than the classical square-root stress singularity. Thus, 62 can be 
either greater or smaller than the conventional inverse square-root 
singularity, depending upon the values of and fiber orientations 

of the composites. Owing to the complex algebraic structure of the 
transcendental characteristic equation, Eq. 20, it is generally not possible 
to predict in explicit form whether 62 > -0.5 or 62 < -0.5 for a delamination 
with crack surfaces in frictional contact without solving the transcendental 
equation numerically. We remark that in the case of a delamlnatlon with crack 
surfaces in frictionless contact between dissimilar anisotropic media, the 
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dominant stress singularity, 6 = -0.5, has also been determined independently 
by using a singular integral equation approach in [10,16]. Furthermore, a 
similar phenomenon of stress singularity 6 > -0.5 or 6 < -0.5 has been 
observed in studying the interface crack between dissimilar isotropic media 
with crack surfaces in frictional contact [28]. 

4.3 Delarainatlon with a Very Small Area of Crack-Tip Closure 

The delamination with open crack surfaces between dissimilar fiber 
composites has been shown mathematically in Section 4.1 to possess 
controversial oscillatory crack-tip stress and displacement fields. This 
abnormality is thought to be artifacts resulting from the method of approach 
by using eigenfunction expansion in the formulation and solution. As first 
pointed out by England [7], Malyshev et al. [6], and later by Wang and Choi 
[10], the region of oscillatory solutions for a delamination with open crack 
surfaces in a nominal tensile field is generally extremely small in comparison 
with the size of the interface crack and this very localized abnormality may 
not be significant in practical terms of linear fracture mechanics. In fact, 
using the partially closed crack model, Wang and Choi [10] have shown that a 
composite delamination in a tensile field has an extremely small crack-tip 
closure with c/a 0(10 A simplified model which disregards the small 

closure (or oscillatory) region and approximates the asymptotic field by an 
inverse square-root stress singularity is, therefore, proposed for this 
situation and shown to provide excellent results comparing with those 
determined by using a partially closed crack model [10]. 

Under certain loading conditions other than pure tension, however, an 
interlaminar crack may also possess a very small area of crack-tip closure, 
depending upon loading modes and material elastic properties of the dissimilar 
composite laminae [16]. A simplified solution for this case can be obtained 
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by taking the limit of the crack-closure length, i.e. , c ->■ 0, in the results 
derived from a partially closed delamination. In particular, the delamlnation 
stress singularity can be taken directly from the partially closed crack model 
(with frictionless crack surface contact) as an inverse square-root one. 
Mathematically, this is equivalent to finding an analytical solution for a 
fully open crack by following the same formulation and procedure for a 
partially closed crack case with infinitesimal closure length, and the 
approximation introduced has the effect of smoothing the oscillatory 
singularity to an inverse square-root singularity for the composite 
delamlnation. Therefore, Interlaminar stress intensity factors and strain 
energy release rates can be defined in a manner consistent with those for a 
homogeneous crack and for the refined model of an interface crack between 
dissimilar isotropic solids introduced by Comninou [14], As will be shown 
later [18], this simplification leads to a very effective and efficient 
approach to the complex problem of delamination with a fully open crack tip or 
with a very small area of crack-tip closure, and provides meaningful 
information on the fundamental mechanics of delamination problems in composite 
laminates under general loading conditions, 

4.4 Logarithmic Stress Singularities 

As mentioned in Section 3 that besides the power-type stress 
singularities given in Sections 4.1, 4,2 and 4.3 for various conditions in the 
crack-tip region, weak logarithmic-type singularities may also appear in the 
homogeneous and particular solutions for the delamination stress field. Since 
the particular solution is related to the remote loading applied to the 
delaminated composite, it has to be considered and constructed for each 
individual case. To study the possible presence of the logarithmic stress 
singularities in a delamination mechanics problem, we consider a symmetric 
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composite laminate subjected to uniform In-plane stretching with 

= Gq for simplicity and without loss of generality. Also, we restrict our 
attention at this point to delamlnatlons located In the following three 
composite systems: 0/-6, 0/0®, and 0/90® graphite-epoxy composites. Based on 

the preceding theoretical developments and the conditions for the presence of 
the logarithmic terms given In [25], we address each Individual case 
separately, 

(1) Delamlnatlon between 0/-0 composites (0^0® and 90®) 

Numerical calculations by using the ply elastic constants given In 
Eq, 28 provide the following: 

N = 12, R = 9, M = 3 (6jj = Integer), (30a) 

N = 12, R = 11, M = 1 (5jj * Integer), (30b) 

Applying both Eq, 30(a) and Eq, 30(b) to the condition 1 = M-(N-R), we obtain 

6 a 

= 0, Thus, logarithmic terms of the form (In Z^) with > 1 do not 

appear in the homogeneous solution for this class of problems. 

Also, carrying out the computations of constructing the left 
*(L) * 

eigenvector C and the loading vector p , we find that, for all three sets 

of ^ at 6 = 0 in this case, Equation 24 is satisfied identically. 

(A^+1) 

Therefore, logarithmic terms of the form (An do not occur in the 

particular solution either. 

(2) Delamination between 0/0® or 0/90® composites (0^0® and 90®) 

In these two cases, following the same procedure and computations as 
discussed in (1) but with minor modifications, we obtain similar results as 
those in the 0/-0 case, i.e., 

N = 12, R = 9, and M =* 3, (6^^ = Integer) (31a) 

N = 12, R = 11, and M = 1, (5^^ ^ Integer) (31b) 
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and, also, = 0. Moreover, Eq. 24 also holds for this problem. Thus, we 

conclude that no logarithmic singularities of any kind would appear in the 

asymptotic solutions for delamination in 0/90° and 0/0° composites; only 

6 

power-type singularities occur in these problems. 

We further remark that, in fact, it has been shown in [25] that no 
logarithmic terms of any kind would occur in the solutions for a general case 
of a delamination located between 0^ and 02 fiber composites with 0j and 0£ 
being any arbitrary fiber orientations. 
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5. SUMMARY AND CONCLUSIONS 

The mechanics of delamination in fiber composite laminates has been 
studied. Formulation of the problem is based on Lekhnitskii 's complex- 
variable stress potentials and basic relationships in laminate elasticity 
theory for anisotropic fiber composites. The eigenfunction expansion method 
used in this study appears to be a suitable approach to determine delamination 
stress singularities and fundamental structures of stress and deformation 
field solutions. Stress singularities for a delamination are found to be 
related to adjacent ply material properties and local traction boundary 
conditions. Numerical results for interlaminar cracks in commonly used 
graphite-epoxy composites with different fiber orientations and crack-tip 
conditions are shown to illustrate the basic nature of stress singularities 
and general solutions for the composite delamination problem. Based on the 
information obtained, the following conclusions may be reached: 

(1) Assuming the delamination is fully open and free from surface 
traction, we find that delamination stress singularities always possess an 
oscillatory form by simultaneous presence of three distinct eigenvalues, 
-V2+iY, -V2-iY. and - V 2 . The oscillatory stress singularities and field 
solutions for composite delamination are physically inadmissible because of 
interpenetration of crack-surface displacements. 

(2) For a delamination with partially closed crack surfaces in 
frictionless contact, the present eigenfunction expansion approach always 
gives an eigenvalue 5 = - V 2 with an algebraic multiplicity of two (l.e., 
double roots), indicating the classical square-root stress singularity is 
recovered in the closed crack case* 
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(3) In the case of a delamination with crack surfaces in frictional 
contact, crack-tip stress singularities depend not only on material elastic 
constants and fiber orientations of adjacent plies but also on frictional 
coefficients fj.^ and along the delamination surface. 

(4) The crack-surface friction may lead to either a stronger or weaker 
stress singularity than the conventional inverse square-root one, depending 
upon fiber orientations of the adjacent plies. Present numerical results, for 
example, show that a weaker stress singularity, l.e., 0 > 6 > ~V2» occurs for 
a delamination between any 0 and -0 fiber composites, but a stronger 
singularity, l.e., -V 2 > ^ > “1» occurs for a delamination between 30® and 
90® composites, if > 0 and ^ 

(5) In the situation that the delamination contains a very small area of 

—6 

crack-surface closure (e.g., c < 10 a), a simplified model with the crack-tip 

stress field having an inverse square -root stress singularity, as determined 
by finding the solution from the limiting case of a partially closed crack 
solution, is suggested and later used for solving the complete boundary value 
problem. 

(6) Examining the multiplicity of eigenvalues and the rank and order of 

the coefficient matrix in the eigenfunction solution, we find that no 

logarithmic stress singularities of any kind would appear in the homogeneous 

and particular solutions for the composite delamination problem; only power- 

5 

type singularities of the form could occur. 

(7) After determining of all the eigenvalues for each individual 
delamination case, general solution structures for composite deformation and 
stress fields can be established immediately. Numerical methods such as the 
singular finite-element technique, which can incorporate exact delamination 
stress singularities in the element formulation, can be easily developed to 
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solve the complete boundary-value problem for delamlnatlons in composite 
laminates with any arbitrary combinations of lamination, geometric, and crack 
variables. One of such methods employing displacement-based singular crack- 
tip elements is given in the associated paper [18], 
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TABLE 1 

DOMINANT STRESS SINGULARITIES* FOR DELAMINATION 
BETWEEN 6/90° GRAPHITE-EPOXY COMPOSITES 


0 



^1 



^2 

^3 

o 

O j 

-0.5 

-h 

0.051101': 

-0.5 

_ 

0.051101: 

-0.5 


(-0.5 

+ 

0.032924:) 

(-0.5 

- 

0.032924:) 

(-0.5) 

15° 

-0.5 

+ 

0.050349: 

-0.5 

— 

0.050349: 

-0.5 


(-0.5 

+ 

0.032484:) 

(-0.5 

- 

0.032484:) 

(-0.5) 

30° 

-0.5 

+ 

0.045138: 

-0.5 

— 

0.045138: 

-0.5 


(-0.5 


0.025764:) 

(-0.5 

- 

0.025764:) 

(-0.5) 

45° 

-0.5 

+ 

0.034504: 

-0.5 

_ 

0.034504: 

-0.5 


(-0.5 

+ 

0.015604:) 

(-0.5 

- 

0.015604:) 

(-0.5) 

60° 

-0.5 

+ 

0.021119: 

-0.5 


0.021119: 

-0.5 


(-0.5 

+ 

0.008067:) 

(-0.5 

- 

0.008067:) 

(-0.5) 

75° 

-0.5 

+ 

0.008899: 

-0.5 

_ 

0.008899: 

-0.5 


(-0.5 

-h 

0.004265:) 

(-0.5 

- 

0.004265:) 

(-0.5) 

90°'*' 

-0 . 5 



-0.5 



-0.5 


(-0.5) 


(-0.5) 



(-0.5) 


*Values in the parentheses are for T300/5280 graphite-epoxy with laminar 
elastic constants given in Eq. 28(a) . 

‘These eigenvalues are for 0°/90° and 90° 790"^ composites in a general loading 

condition. In the cases of 0°/90° and 90°/90° composites under uniform 

stretching £_=£_, 6o = -0.5 does not appear because of = x =0 and 
u ' ^ xz y z 

being decoupled from other stress components. 
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TABLE 2 

DOMINANT STRESS SINGULARITIES* FOR DELAMINATION 
BETWEEN 30° /e GRAPHITE-EPOXY COMPOSITES 


0 



^2 

s 

0° 

-0.5 + 0. 01245 
(-0.5 + 0.009456i) 

-0.5 - 
(-0.5 - 

0.01245U 

0.009456i) 

-0.5 

(-0.5) 

15° 

-0.5 + 0.010491i 
(-0.5 + 0.007920^) 

-0.5 - 
(-0.5 - 

0.010491i 

0.007920i) 

-0.5 

(-0.5) 

30° 

-0.5 

(-0.5) 

-0.5 

(-0.5) 


-0.5 

(-0.5) 

45° 

-0.5 + 0.015968i 
(-0.5 + 0.011608^) 

-0.5 - 
(-0.5 - 

0.015968i 

0.011608i) 

-0.5 

(-0.5) 

o^ 

o 

o 

-0.5 + 0.030943i 
(-0.5 + 0.020441i) 

-0.5 - 
(-0.5 - 

0.030943i 

0.020441i) 

-0.5 

(-0.5) 

75° 

-0.5 + 0.041030i 
(-0.5 + 0.024664i) 

-0.5 - 
(-0.5 - 

0.041030i 

0.024664i) 

-0.5 

(-0.5) 

90° 

-0.5 + 0.045138i 
(-0.5 + 0.024764i) 

-0.5 - 
(-0.5 - 

0.045138i 

0.025764£) 

-0.5 

(-0.5) 


Values in the parentheses are for T300/5280 graphite-epoxy with laminar 
elastic properties given in Eq. 28(a) • 
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TABLE 3 

DOMINANT STRESS SINGULARITIES* FOR DELAMINATION 
WITH CRACK-TIP CLOSURE IN 0/-0 GRAPHITE -EPOXY COMPOSITES 


15“ 30“ 45“ 60“ 75“ 


0.0 

-0.5 

- 0.5 

- 0.5 

-0.5 

- 0.5 


-0.5 

- 0.5 

- 0.5 

-0.5 

- 0.5 

0.2 

-0.5 

- 0.5 

- 0.5 

-0.5 

- 0.5 


-0.4994 

-0.4962 

-0.4940 

-0.4941 

-0.4964 

0.4 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.4977 

-0.4923 

-0.4881 

-0.4881 

-0.4928 

0.6 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.4965 

-0.4884 

-0.4821 

-0.4822 

-0.4892 

0.8 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.4954 

-0.4846 

-0.4762 

-0.4763 

-0.4855 

* 6 ]^ and 

62 are found to 

be independent of the value 
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TAB LI' 4 

DOMINANT STRESS SINGULARITIES FOR DELAMINATION 
WITH CRACK-TIP CLOSURE IN 30°/90° GRAPHITE/EPOXY COMPOSITES 


- 

• r(? \ 

0.0 

0.1 

0.3 

0.5 

0.7 

0.0 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.5 

-0.5010 

-0.5031 

-0.5051 

-0.5072 

0.1 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.5036 

-0.5046 

-0.5067 

-0.5087 

-0.5108 

0.3 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.5108 

-0.5118 

-0.5138 

-0.5159 

-0.5179 

0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.5179 

-0.5189 

-0.5210 

-0.5230 

-0.5251 

0.7 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 


-0.5251 

-0.5261 

-0.5281 

-0.5302 

-0.5322 
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7, LIST OF FIGURE CAPTIONS 


Fig. 1 Delamlnations In a [01/02/63/. . ./03/02/01I Fiber-Reinforced Composite 
Laminate . 

Fig. 2 Coordinates and Geometry of a Delamination with Open Crack Surfaces 
between 0^^ and ®m+l Plies. 

Coordinates and Geometry of a Delamination with Finite Length of 
Crack-Surface Closure between 0j^ and 6^+1 Plies. 


Fig. 3 





Fig, 3 Coordinates and Geometry of a Delamination with Finite Length of 
Crack-Surface Closure between 9^ and 0„+l Plies. 


O 
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APPENDIX 1 


Expressions for coefficients in Eqs. 15 and 16 are as follows: 


^Ik “ coS(f))2, 


H 2 j^= -nj^(yj^sln<|» + cos(|)). 


= -(Pj^sin(|) + cos(())(Uj^cos(j) - sin(J)), = (pj^cos<|) - sin<|))2, (A-1) 


= t|j.(Mj^cos<|) - sin(|)). 


^6k " Pk*^®®*** 


Hyk = -Pj^^sin(|) + qj^cos(|). 


®8k “ *^k» 


where \i^ and are defined in Eqs* 11, 12 and 13, and pj^, q|^ and are 
complex constants related to laminar elastic constants by 

\ “ ^ll^k ^12 “ ^14\ ^15W “ ^le^^k’ 


\ “ ^IZ\ ^ hlK ~ ^ 24 \/\ ^ 25 \ " ^ 26 ’ 


(A-2) 


*^k ~ ^lA’^k ~ ^4aV\ ^45\ “ ^46* 
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